Modified Fermi-sphere, pairing gap and critical temperature for the BCS-BEC 
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We investigate the phase diagram of two-component fermions in the BCS-BEC crossover. Using 
functional renormalization group equations we calculate the effect of quantum fluctuations on the 
fermionic self-energy parametrized by a wavefunction renormalization, an effective Fermi radius and 
the gap. This allows us to follow the modifications of the Fermi surface and the dispersion relation 
for fermionic excitations throughout the whole crossover region. We also determine the critical 
temperature of the second order phase transition to superfluidity. Our results are in agreement 
with BCS theory including Gorkov's correction for small negative scattering length a and with 
an interacting Bose gas for small positive a. At the unitarity point the result for the gap at 
zero temperature agrees well with Quantum-Monte-Carlo simulations while the critical temperature 
differs. 

PACS numbers: 03.75.Ss; 05.30.Fk 



I. INTRODUCTION 

In the vicinity of a Feshbach resonance cold fermionic 
atom gases show a continuous crossover from BCS to 
BEC condensation [l|. This qualitative picture has been 
confirmed in recent experiments As the system can 
serve as a testing ground for non-perturbative methods 
a quantitatively precise theoretical understanding of this 
phenomenon is the aim of various theoretical approaches. 
At and near the resonance quantum Monte-Carlo (QMC) 
results have been obtained The complete phase di- 

agram has been accessed by functional field-theoretical 
techniques, in particular by e-expansion [8;ll[, 
expansion [l^ [l3| , i- matrix approaches Il5|. Dyson- 
Schwinger equations [Hi, J7,], 2-Particle Irreducible meth- 
ods (isj , and nonp erturbative renormalization-group flow 
equations [19l424| . 

The techniques and results from renormalization-group 
flow equations for cold atoms have been gathered and 
reviewed in psj . There, a basic truncation is investi- 
gated and extended to include important quantitative 
effects, such as particle-hole fluctuations. In the present 
work we aim at a systematic completion of the trun- 
cation scheme and consider further important, however 
subleading terms in the flowing action. In the fermionic 
sector we include a renormalization term adding to the 
chemical potential and the fermionic wave function renor- 
malization. This accounts for changes of the Fermi sur- 
face due to fluctuation effects. Similar parametrizations 
of the fermionic self-energy have been studied in refs. 
[Hill Hi- In ref. f??] the spectral function of ultracold 
fermions in the BCS-BEC crossover has been calculated 
from a ladder approximation to the Luttinger-Ward func- 
tional. 

As an additional parameter, we also include an atom- 
dimer interaction term. These extensions of the trunca- 
tion improve our numerical results in the perturbative 



domains of the crossover to match the results from well- 
known other fleld-theoretical methods. We obtain con- 
siderable quantitative precision on the BCS and BEC 
sides of the resonance. Also in the strongly interacting 
regime of the Feshbach resonance we obtain a good agree- 
ment with data from QMC simulations for the gap and 
the chemical potential at zero temperature. The effect 
of the additional couplings on the critical temperature 
at the unitarity point, however, is rather modest and no 
agreement with the widely accepted QMC simulations is 
found. 

For a microscopic model of the BCS-BEC crossover we 
employ a two-component Grassmann field ^ = (V'1,'02), 
describing fermions in two hyperflne states. As bosonic 
degrees of freedom we introduce a complex scalar field 
(j), which can describe molecules, Cooper pairs or simply 
an auxiliary field in different regimes of the crossover. 
The resulting two-channel model can describe both nar- 
row and broad Feshbach resonances in a unified setting. 
The microscopic action at the ultraviolet scale A reads 
explicitely 



l/T 



dr / d-^x 



-/ia((?!)*?AiV'2 + h.c.)^ , 



(1) 



where we choose nonrelativistic natural units with h = 
ks ~ 2Af = 1, with M the mass of the atoms, and A 
denotes the Laplace operator. We use the Matsubara for- 
malism to describe the system in thermal equilibrium, so 
the fields depend on the position variable x and on the 
imaginary time variable r, which parameterizes a torus 
with circumference l/T. The variable /i is the chemical 
potential. We introduce a Yukawa or Feshbach interac- 
tion h which couples the fermionic to the bosonic fields 
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and directly relates to the width of the Feshbach reso- 
nance. The parameter v depends on the magnetic field 
and accounts for the detuning from the Feshbach reso- 
nance. The microscopic values /ia, i^a can be fixed by 
the properties of two bod y sc attering in vacuum. This is 
discussed in detail in (20l. l25l|. 

Here we will discuss the limit of broad Feshbach res- 
onances, which corresponds to /ia — > oo. Then the mi- 
croscopic interaction between fermions becomes point- 
like, with strength —h\/vf^. The system shows a far 
reaching universality at a broad Feshbach resonance such 
that macroscopic quantities are independent of the mi- 
croscopic details. In our case the relevant parameter is 
the two-body scattering length a or, at finite density, the 
concentration c — akp, where the formal Fermi momen- 
tum is related to the density by kp — (Stt^ti)^/^. At 
nonzero temperature we have an additional parameter 
given by T jTp, where Tp = Ep = kp is the Fermi tem- 
perature. 

The paper is organized as follows: In Sect. [ll] we in- 
troduce the exact functional flow equation and specify 
the truncation as well as our regulator scheme. Sec- 
tion mil discusses the treatment of particle-hole fluctu- 
ations by a scale-dependent Hubbard-Stratonovic trans- 
formation using a recently derived exact flow equation. 
We present our results for the fermionic self-energy, the 
single-particle gap and the critical temperature through- 
out the crossover in Sect. IIVI and finally draw some con- 
clusions in Sect. IVl 



II. FLOW EQUATION & TRUNCATION 

The functional renormalization group connects the 
microphysics, cf. equation ([T]), to macrophysics, i.e. 
to observable thermodynamics, by means of a non- 
perturbative flow equation. The scale dependence of the 
flowing or average action is given by the exact functional 
fiow equation j28|] 



dkRk 



(2) 



Here, the STr operation involves an integration over mo- 
menta as well as a summation over internal indices, with 
an appropriate minus sign for the fermions. The infrared 
cutoff function Rk{q) has the properties Rk{q) -> oo for 



A, Rkiq) 



k^ for 



l9l 



and Rk{q) for 



ill — T' 0. We employ the second functional derivative of 

Tk 



rr^M^^- (3) 



Equation ^ is an exact differential equation for a 
functional. This corresponds to a system of infinitely 



many coupled differential equations for running cou- 
plings. Finding exact solutions to equation ^ for non- 
trivial theories is possible in practice, but one can use 
truncations in the space of possible functionals to find 
approximate solutions. Such approximations do not have 
to rely on the existence of a small expansion parameter 
like the interaction strength and they are therefore of a 
non-perturbative nature. For reviews of the functional 
renormalization group method see [H, H^l • 

In this work, we solve the flow equation ^ by us- 
ing a truncation which is a completion of the truncation 
scheme that was discussed in detail in [2^, [2^. More 
explicitly. 



l/T 



(4) 



The effective potential U (p, /i) contains no derivatives 
and is a function of p = (f>*(f> and /i. Besides the couplings 
parameterizing U (see below) our truncation contains six 



further fc-dependent ("running") couplings Z. 



A, 



Z^, h, and A^^. In addition to the parameters used in our 
previous investigations [23, 25], our truncation includes 
now the fc-dependent parameters and Z^ by which 
we take into account fluctuation effects on the self-energy 
of the fermionic quasiparticles. The inclusion of the cou- 
pling constant A^^ closes the truncation on the level of 
interaction terms quartic in the fields (see also [3l|). The 
coupling Aw,,/, plays an important role for the three-body 
physics [32, ^53] and is expected to lead to quantitative 
modifications for the many-body problem. In contrast to 
[3^ [33I we do not resolve the momentum dependence of 
\tj,^ in this work. 

Another quantitative contribution to the many-body 
physics of the BCS-BEC crossover is given by particle- 
hole fiuctuations. We will include them by the technique 
of scale dependent bosonization. A brief discussion of 
this is given in the next section. For a more detailed 



discussion we refer to a previous paper 23 1 



In terms of renormalized fields 

yl/2 



P = A^p, 



h = 



ip = Z^ -0, renormalized couplings Z^ — Z^ 
hl{A^J'^ Z^), A^^ = \^,p/{A^Z^), m'^ ^ ffi'^/Z^ and ef- 
fective potential U{p,p) — U{i5,p), Eq. ^ reads 



+ r{Z^dr-^A)^ + U{p,p) 



- /i(0*V^iV'2+0V'2>D+A0^0*0V'V|- (5) 

For the effective potential, we use an expansion around 
the fc-dependent location of the minimum pQ{k) and the 
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fc-independent value of the chemical potential fiQ that 
corresponds to the physical particle number density n. 
We determine po(fc) and by 



{dpU){po{k),^io) 







for all k 
at fc = 0. 



(6) 



More explicitly, we employ a truncation for U (p, /x) of the 
form 



U{p,p.) U{po,po) -Ukip- p,o) 

+ {m^ + a{p, ~ Mo))(p - Po) 



(7) 



This expansion is discussed in detail in refs. [3J, |35|. In 
the symmetric or normal gas regime, we have pq = 0, 
while in the regime with spontaneous symmetry break- 
ing, we have = 0. The atom density n = —dU /dpi 
corresponds to rifc in the limit fc 0. The system is su- 
perfluid if pQ{k — !• 0) > 0. For po > the term ~ A^^ 
leads to a further modification of the Fermi surface, be- 
sides the gap \fWpQ. 

In total, we have eight running couplings m'^{k), 
m'^{k), A(fc), a{k), Uk, Z^{k), h{k) and \^^{k). (In 
the phase with spontaneous symmetry breaking m? is 
replaced by po-) In addition, we need the anomalous 
dimensions rj^ — — fcc^fclnA^ and ry^ = — fc9fclnZ^. We 
project the flow of the average action on the flow 
of these couplings by taking appropriate (functional) 
derivatives on both sides of Eq. ^ . Wc thereby obtain a 
set of coupled nonlinear differential equations which can 
be solved numerically. 

At the microscopic scale fc = A the initial values of 
our couplings are determined from Eq. ([1]). This gives 
m|(A) = -p, m^{A) = z^a - 2p, po,A = 0, A(A) = 0, 
Z^{A) = 1, h{A) = Ha, a(A) = -2 and riA = 3T:^pe{p). 
The initial values va and can be connected to the 
two particle scattering in vacuum close to a Feshbach 
resonance. For this purpose one follows the flow of rn^(fc) 
and h(k) in vacuum, i.e. p, = T — n — and extracts the 
renormalized parameters = m?(k = 0), h = h{k = 0). 
These are connected to the scattering length by 



h^{k = 0) 
8TTm^{k = 0)' 



(8) 



For details we refer to ref. [25]. 

The infrared cutoff we use is purely space-like and 
reads in terms of the bare fields 



Pk 



-iJj\p)Z^ [sign(p*^ ~Pk)k^ 

{f-~Pk))]e{k^-\f-~Pk\)i^{p) 

4>*{p)A^ [e - f/2\ e {e - f/2) m 

— TO?, — X^tpPO- 



For the fermions it regularizes fluctuations around the 
running Fermi surface, while for the bosons fluctuations 
with small momenta are suppressed. The choice of ASV 
in Eq. ^ is an optimized choice in the sense of [13, HH, 



III. SCALE DEPENDENT BOSONISATION 

In the microscopic model in Eq. ([T]) a Hubbard- 
Stratonovich transformation relates the model for 
fermions with a local interaction to our Yukawa-like the- 
ory. The fermion-fermion interaction has been replaced 
by a tree-level process that involves the composite bo- 
son (j). In the limit of vanishing temperature and density 
(few-body physics) this leads to a consistent treatment 
since an explicit fermion-fermion interaction term in the 
effective action vanishes on all scales. At nonzero density 
the situation is different, however. Particle-hole fluctua- 
tions generate a term of the form 



l/T 



drdPx '01 "01 ■02 "02 



(10) 



in the effective action, although it vanishes in the micro- 
scopic action, cf. Fig. [T] In principle one could simply 
take A^ as an additional coupling into account. How- 
ever, it is much more elegant to use a scale-dependent 
Hubbard- Stratonovich transformation [ssl . which ab- 
sorbes A^ into the Yukawa-type interaction with the 
bosons at every scale fc. By construction, there is then 
no self interaction between the fermionic quasiparticles. 
The general procedure of "partial bosonization" is dis- 
cussed in detail in [23l |. In the present paper we employ 
the exact flow equation derived in [39']. Compared to 
the scheme used in [23i] the flow equations are modified 
slightly. Following the steps in [23, [s^ we find for the 
renormalized coupling to^ in the symmetric regime an 
additional term reflecting the absorption of A^ into the 
fermionic interaction induced by boson exchange. 



HS 



/l2 



■dtX 



(11) 



Here c^jm^jj^g and 9(A^i,|jjg denote the flow equations 
when the Hubbard-Stratonovich transformation is kept 
fixed. The flow equation of all other couplings are the 
same as with flxed Hubbard-Stratonovich transforma- 
tion. In the regime with spontaneous symmetry breaking 
we use 



dth 

dtPo 
dtX 



dtX\ 



HS 



(12) 



(9) 
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A. Fermi sphere 



dtX^ 
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FIG. 1: Box diagram for the flow of the four-fermion interac- 
tion. 
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FIG. 2: Flow of the inverse fermionic wave function renor- 
malization at T = 0. The solid line is at the unitarity 

point (c~^ = 0). The dashed line corresponds to the BCS 



side at c 



at c 



= 1. 



-1 and the short dashed hne to the BEG side 



IV. RESULTS 



The extension of the truncation by including the run- 



ning couplings Z. 



and Aa,/, and the different bosoni- 



sation scheme as compared to [23| help to improve the nu- 
merical precision of our results considerably all over the 
crossover, however especially in the limits of weak inter- 
actions. The effect of including the A^^-vertex can most 
prominently be observed in the result for the vacuum 
dimer-dimer scattering length ajvf/a. In this truncation 
we find ajv/Za = 0.59, which is in very good agreement 
with the well-known result from solving the Schrodinger 
equation, auja. — 0.60 40]. Our result is somewhat sur- 
prising since no momentum dependence of A^^ has been 
taken into account while this turns out to be quite im- 
portant for the atom-dimer scattering length [32|. 

A study of the flow of the wave function renormaliza- 
tion Zy, shows, that its effect is strongest in the strongly 
interacting regime, when |c~^| < 1, cf. Fig. [2] Away 
from the unitarity point has only a weak scale depen- 
dence and stays close to its initial value Z^(k = A) 1, 
see Fig. [31 



parametrize the effect of fluctuations on the fermionic 
self-energy as well as the coupling A^^ have an interest- 
ing effect on the dispersion relation of the fermions. In 
the regime with spontaneously broken C/(l) symmetry, 
the renormalized propagator of the fermionic field reads 
after analytic continuation to real frequencies w 



G- 



-h4)Qe, 



-uj - + m'^ + A^^po) 



h<j)oe 



(13) 



The dispersion relation follows from det ^ = as 



(14) 



where A = /ly^ is the gap and rp = y ""^^ ~ X^tpPo is 
the effective radius of the Fermi sphere. For small nega- 
tive scattering length on the BCS side of the crossover the 
renormalization effects on ri? remain small. It takes the 
classical value rp = ^ = kp where kp = (Stt^m)^/'^. 
For large negative scattering length the Fermi sphere 
gets smaller and rp finally vanishes at a point with 
{akp)~^ K, 0.6. For small and positive scattering length 
(BEC side), the Fermi surface has disappeared and the 
fermions are gapped even for A — > by a positive value 
of TO^ + A^^po- We plot our result for the Fermi radius 
in units of the Fermi momentum kp = (Stt^tt.)^/'^ in Fig. 
|4]and the dispersion relation in Fig. [5l 

In order to show explicitely the effect of the atom- 
dimer vertex A^^ on the radius of the Fermi sphere we 
also plot rp as A^^ is omitted in the truncation (Fig. UJ 
dotted line). We find that this changes the estimate for 
rp/kp quantitatively for crossover parameters (akp)~^ > 
— 1. We conlude that the atom-dimer vertex affects not 
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FIG. 3: Inverse fermionic wave function renormalization 1/Z^ 
at the macroscopic scale fc = as a function of the crossover 
parameter [akp)'^ ■ 
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FIG. 4; Effective Fermi radius rp/kp as a function of the 
crossover parameter (akF)~^ for vanishing temperature (solid 
line). We find that the Fermi sphere vanishes approximately 
at the point with {akF)~^ ~ 0.6. For comparison we also 
give the effective Fermi radius in an approximation, where 
we omitted the contribution of the atom-dimer vertex A^^ 
(dotted line). 



uj/E_ 




FIG. 5: Positive branch of the dispersion relation u}{q) in 
units of the Fermi energy as given in Eq. (|14|) for — —1 
(dashed), c^^ = (solid) and c^^ = 1 (short dashed). The 
gap and the Fermi radius are evaluated at the macroscopic 
scale fc = 0. 




FIG. 6: Gap in units of the Fermi energy A/Ep as a func- 
tion of (akp)^^ (solid line). For comparison, we also plot 
the result found by Gorkov and Melik-Bakhudarov (left da- 
hed) and extrapolate this to the unitarity point {akF)~^ = 0, 
where Aqmb/Ef — 0.49. 



lated from perturbation theory. The result first found 
by Gorkov and Melik-Bakhudarov is given by 



A/Ef 



2X 7/3 



(15) 



In our approach we can extend this calculation to the 
strongly interacting regime and to the BEG side of the 
crossover. Our result for the gap in units of the Fermi 
energy, A/Ep, is shown in Fig. [51 We find perfect agree- 
ment with Eq. (|15l) on the BGS side. At the unitarity 
point, (afci?)~^ ~ 0, we obtain A/Ep — 0.46. Another 
benchmark for the comparison of different methods is 
the chemical potential in units of the Fermi energy at 
(akp)-'^ = and T = 0. We find (m/Ef = 0.51 and com- 
pare the results of different approaches to this quantity 
as well as A/Ep in Tab. HI 



C. Critical temperature 

We next turn to the result for the critical temperature 
throughout the whole crossover which is shown in Fig. [T] 



only the physics on the BEG side as reflected for example 
in the ratio / o-, but also has an effect close to unitarity. 



B. Single-particle gap at T = 

To reduce the numerical effort we have omitted the 
effect of particle-hole fluctuations for the study of the 
effects discussed so far. In the following all results are 
given inluding particle-hole fluctuations, however. 

On the BGS side for small negative values of akp the 
single particle gap at zero temperature can be calcu- 



TABLE 1: Results for the single-particle gap and the chemi- 
cal potential at T = and at the unitarity point by various 



authors. 

fi/EF A/Ef 

Carlson et al. [3] 0.43 0.54 

Perali et al. [15| 0.46 0.53 

Haussmann et al. fl8] 0.36 0.46 

Diehl et al. [20] 0.55 0.60 

Bartosch et al. [24] 0.32 0.61 

present work 0.51 0.46 
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On the BCS side (c ^ < 0) we observe a perfect matching 
with BCS theory including the correction by Gorkov and 
Mehk-Barkhudarov El 



^ = 0.28e"/(2c) 
Tf 



(16) 



which is given by the short dashed hne in Fig. [71 De- 
viations from this pertubative regime appear only quite 
close to the regime of strong interactions —> 0. 

On the BEC-side our result approaches the critical 
temperature of a free Bose gas. This value is approached 
in the form [H 



Tc — Tc BEC 0.M 
~ K 



T, 



c,BEC 



a (67r2)i/3' 



(17) 



as shown by the long dashed line in Fig. [T] Here, qm 
is the scattering length between the molecules. For the 
ratio aM/a we use our result aM /d = 0.59 obtained from 
solving the flow equations in vacuum. This has to be 
compared to the result obtained from solving the corre- 
sponding Schroedinger equation which gives aujo- = 0.6 
[40] ■ For the coeffcient k determining the shift in Tc com- 
pared to the free Bose gas we find k = 1.39. In [H, |43| 
the result for an interacting BEG is given with k = 1.31, 
which is in reasonable agreement with our result. Other 
characteristic quantities are the maximum of the ratio 
(Tc/rp)inax ~ 0.31 and the location of the maximum 
(a^F)-L ~ 0.40. 

In the unitary regime (c ^ « 0) our extensions of the 
truncation have a considerable quantitative effect. We 
can give an improved estimate for the critical tempera- 
ture at the resonance (c~^ = 0) where we find Tc/Tp = 




FIG. 7: Critical temperature Tc/Tp in units of the Fermi 
temperature as a function of the crossover parameter (akp)'^ 
(solid line). For comparison, we plot the BCS result with 
Gorkov's correction (left short dashed). On the BEC side we 
show the critical temperature for a gas of interacting bosonic 
molecules according to [i^ - li^ (dashed line on the right) . The 
three red dots close to and at unitarity show the QMC results 
by Burovski et al. [3]. 



0.248 and a chemical potential ^c/Tp = 0.55. This com- 
pares to the previous FRG estimate T^jTp = 0.264 and 
IJ-c/Tp — 0.67 in ref. ,23]. Results from Quantum Monte 
Carlo simulations are T^/Tp — 0.15 and Hc/Tp = 0.49 
in '4], Tc/Tp < 0.15 and ^ic/Tp = 0.43 in Q and 
Tc/Tp = 0.245 in 0. 



V. CONCLUSIONS 

We have studied the flow equation approach to the 
BGS-BEG crossover in a rather detailed truncation. In 
addition to our previous treatments we have included 
contributions to the fermionic self-energy, a further in- 
teraction between fermions and bosons and improved our 
treatment of particle- hole fluctuations. Our results are in 
very good agreement with analytical treatments in vari- 
ous regimes. 

For small negative values of the concentration parame- 
ter akp our findings for the gap at vanishing temperature 
A/Ep and the critical temperature Tc/Tp confirm the 
predictions of Gorkov and Melik-Barkhudarov [4l|. On 
the other side, in the BEC regime with small positive val- 
ues oi akp, our results are in agreement with the expecta- 
tions for a weakly interacting Bose gas where the shift in 
the critical temperature is linear in the dimer-dimer scat- 
tering length um- Also our finding for the ratio am /a is in 
goo d agreement with the quantum-mechanical treatment 
[40| . This quantitative accuracy is remarkable in view of 
the fact that we have started with a purely fermionic mi- 
croscopic theory, without propagating bosonic degrees of 
freedom or bosonic interactions. 

We present results for the single particle gap and the 
modifications of the Fermi surface at T = over the 
whole range of the crossover. This allows us to compute 
the dispersion relation for the single particle or fermionic 
excitations. In the strongly interacting regime where the 
scattering length diverges, no analytical treatments are 
available. Our results for the gap A/Ep and the chemical 
potential fJ./Ep at zero temperature are in reasonable 
agreement with Monte-Carlo simulations. This holds also 
for the ratio fic/Ep at the critical temperature. The 
critical temperature Tc/Tp itself is found to be larger 
than the widely accepted Monte-Carlo result, however. 

In future studies our approximations might be im- 
proved mainly at two points. One is the frequency- and 
momentum dependence of the boson propagator. In the 
strongly interacting regime this could be rather involved, 
developing structures beyond our current approximation. 
A more detailed resolution might lead to modifications 
in the contributions from bosonic fluctuations to vari- 
ous flow equations. Another point concerns structures 
in the fermion-fermion interaction that go beyond a di- 
atom bound-state exchange process. Close to the unitar- 
ity point, other contributions might arise, for example in 
form of a ferromagnetic channel. While further quanti- 
tative modifications in the unitarity regime are conceiv- 
able, the present status of approximations already allows 
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for a coherent description of the BCS-BEC crossover for under FOR 723, and GK1523/1. 
all values of the scattering length, temperature and den- 
sity by one simple method and microscopic model. This 
includes the critical behaviour of a second order phase 
transition as well as the temperature zero range. 
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